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Abstract. In this paper we consider the model of n non-intersecting squared 
Bessel processes with parameter a, in the confluent case where all particles 
start, at time t = 0, at the same positive value x = a, remain positive, and 
end, at time T = t, at the position x = 0. The positions of the paths have a 
limiting mean density as n — > oo which is characterized by a vector equilib- 
rium problem. We show how to obtain this equilibrium problem from different 
considerations involving the recurrence relations for multiple orthogonal poly- 
nomials associated with the modified Bessel functions. 

We also extend the situation by rescaling the parameter a, letting it increase 
proportionally to n as n increases. In this case we also analyze the recurrence 
relation and obtain a vector equilibrium problem for it. 



1. Introduction and statement of results 

1.1. Introduction. This paper deals with the model, studied in 1 2J , of n non- 
intersecting squared Bessel paths in the confluent limit where all paths start, at 
time t = 0, at the same positive value x = a, remain positive, and end at a later 
time t = T, at the position x = 0. 

The squared Bessel process is a diffusion process on [0, oo), depending on a 
parameter a > — 1, whose transition probability density is given by 



P f a (0,y) = r J , -e~*, y>0. 

Here, I a denotes the modified Bessel function of the first kind of order a, 

I ( z ) = V ( ' ' 

la W -2s kWtk + a + lY 

If we let n — > oo, and perform an appropriate time scaling t h-» t/(2n), T t— > 
l/(2n), the paths fill out a region in the ta-plane as shown in the left figure in 
Figure [TJ The figure shows a numerical simulation of 50 non-intersecting squared 
Bessel paths and the boundary of the region filled by the paths. The right part 
of Figure [T] shows a similar picture where in addition a increases proportionally 
with ro. Since a is a measure for the repulsion from all paths in the right part of 
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Figure 1. Numerical simulation of n — 50 non-intersecting 
squared Bessel paths. The left hgure shows n paths with fixed 
a, and the right figure shows n paths with a increasing propor- 
tionally with n. 



Figure [T] stay at a positive distance from until the final time t = 1. In contrast, 
in the left part of Figure [TJ where a remains fixed, the smallest paths arrive at the 
hard edge (the wall) at x — already at a critical time t = t* — a/(l + a) < 1, see 

The model with fixed a was studied in [12]. One of the results was that for 
every t G (0, 1), the positions of the paths at time t have a limiting mean density 
as n — > oo, that is characterized by the following vector equilibrium problem (see 
Theorem 2.4 and the appendix of |12|). The vector equilibrium problem asks to 
minimize the functional 



(1.1) 



log 



1 



\x-y\ 



du 1 {x)dv 1 (y) 



log 



1 



\x - y\ 



dv 2 {x)dv 2 (y) 



log- 



1 



-dv x (x)dv 2 {y) 



\x - y\ 

over all vectors of measures (y\, v 2 ) such that 



1\fax 



t(l-t) t 



dv\{x), 



(1.2) supp(^) C [0, oo), fdv 1 = l, supp(z/ 2 ) C (-oo,0], Jdv 2 = l/2, 



and 
(1.3) 



v 2 < a, 



where a is the measure on (— oo,0] with density 



(1.4) 



da{x) 



-1/2 



X G 



-oo,0]. 



dx nt 

There is a unique minimizer (y\,v 2 \ and the part v\ has a density which is the 
limiting mean density of the positions of the paths at time t. The equilibrium 
problem depends on the parameters a > and < t < 1 that appear in the energy 
functional fll.ip as well as in the upper constraint (|1.4[) . 

The above equilibrium problem is somewhat unusual since it involves a second 
measure v 2 with a constraint a. It turns out (see [12] ) that the constraint is not 
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active if t < t* = a/(a + 1), that is before the critical time where the smallest paths 
come to the hard edge. The constraint is active for t > t* . 

It is the aim of this paper to give more insight into the nature of the vector 
equilibrium problem with constraint. We show how it arises from different con- 
siderations involving the recurrence relations for multiple orthogonal polynomials 
associated with modified Bessel functions, see the next section. 

We note that a similar vector equilibrium problem with constraint and external 
field was also found in the context of the Hermitian two-matrix model 

— exp (-nTr(V(Afi) + W(M 2 ) - tMxM 2 )) dM x dM 2 

where W{M) — jM 4 , see [7J. The vector equilibrium problem now involves three 
measures, which arc located on M, iM. and R, respectively, with a constraint acting 
on the second measure. The energy functional is similar to (jTTTJ). The polynomials 
that are connected to this model satisfy a five term recurrence relation. The ideas 
that are developed in this paper can be extended to analyze this recurrence relation 
and it is possible to obtain the vector equilibrium problem from it. This will be 
reported elsewhere. 



1.2. Multiple orthogonal polynomials. In [TJ] it was shown that the positions 
of the paths at any time t € (0, T) constitute a multiple orthogonal polynomial 
(MOP) ensemble associated with the two weight functions 



(1.5) wi{x)=x^e ^t 1 "*) I a 



t 



W 2 (x) = X 



t 



defined for x G (0, +oo). It implies in particular that the 'average characteristic 
polynomial' 



(1.6) 



B n (x)=E 



Y[(x- Xj {t)) 

i=i 



where x\(t) < X2{t) < ■ ■ ■ < x n (t) are the positions of the n paths at time t, satisfies 
the orthogonality relations (assume n is even) 



(1.7) 



B n (x)x k w 3 (x)dx = 0, for k = 0, 1, . . .,n/2 - 1, j = 1, 2. 



The polynomial (| 1 . 6[) is characterized by a 3 x 3 matrix valued Riemann-Hilbert 
problem and the correlation kernel for the MOP ensemble can be written directly in 
terms of the solution of the Riemann-Hilbert problem. See also [TT] for the general 
notion of a MOP ensemble. 

The multiple orthogonal polynomials B n for the weights (|1.5[) were studied in 
detail by Coussement and Van Assche [1], [5]. They obtained a third order differ- 
ential equation, which was used in [12j . In addition they also gave an explicit four 
term recurrence relation 



(1.8) xB k (x) = B k+1 {x) + b k B k (x) + CkBk-^x) + d k B k ^ 2 {x) 
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with recurrence coefficients 

_ a(T - t) 2 2t(T - t) 

°k — rj, 2 T 



(2k + a + l), 



4 *(fc-l) 



(1-9) c fc = V y3 fc + y T2 k{k + a), 

Aat 2 (T -t) A , „ 

d* = ^^fc(fc-i). 

After the rescaling i i— > t/(2n), T i— > l/(2n), we obtain a doubly indexed sequence 
i?fc,n of polynomials satisfying for each n, the recurrence 

(1.10) xB k , n (x) = B k+ i >n (x) + b kin B k , n (x) + c k , n Bk-i,n{x) + d k , n B k ^ 2 ,n{x) 
with recurrence coefficients 

b k ,n = 0(1 - * 2 + t 1 - t) , 

n 

(1.11) c fen = 2ot(l - i) 3 - +i 2 (l - <) 

71 

d fc ,„ = at 2 (l-t)' 

It is the aim of this paper to obtain the vector equilibrium problem from this recur- 
rence relation. We will show that the measure v\ is the limiting zero distribution 
of the diagonal polynomial B n n as n — + oo. 

We also consider the following extension of the situation studied in [12] . We 
rescale the parameter a by letting it increase proportionally to n as n increases. 
The right figure in Figure [T] shows a numerical simulation of 50 non-intersecting 
paths for this case. We observe that now the paths stay away from the hard edge 
x = 0. In this case we are also able to analyze the recurrence relation and obtain 
a vector equilibrium problem from it. 

1.3. Polynomials satisfying an m-term recurrence relation. The recurrence 
coefficients ([1.111) are such that whenever n — > oo, k — > oo so that k/n — > s exists, 
there is a limit 

b k , n -> b(s), c k ,n -* c(s), d k>n -> d(s) 

for certain functions b(s), c(s) and d(s). 

We will study this situation in the general context of polynomials P k>n depending 
on two parameters, satisfying for each n, an m-term recurrence relation with varying 
coefficients 
(1.12) 

xP k , n (x) = Pk+l, n (x) + bf n P k , n {x) + b^Pk-l,n( x ) + ' ' ' + b k. n % 'flt+2-m,n(a;), 

with Po = 1, P-\ = 0, . . . , P- m +2 = 0, where the recurrence coefficients have 
scaling limits 

lim b k j) n = b {j) (s), j = 0, . . . , m - 2. 

A:/n — >s ' 

for certain functions . . . , b^ m ~ 2 \ The notation limw n _, s that we use here and 
also later in the paper, means that both k, n — > oo with k/n — > s. 

We associate with the functions b^\ j = 0, . . . , m — 2, a family of functions 

(1.13) A s (z) = z + 6 (0) (s) + b {1 \s)z- 1 + ■■■ + b {m - 2 \s)z- m+2 
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and the sequence of Toeplitz matrices (T n ) n where T n = T n (A s ) is the nxn Toeplitz 
matrix with symbol A s , defined by 

fl, Xk = j+1, 

(1.14) (T n {A s )) jk = I b®(a), if A: = j — i, i = 0,...,m-l, 

I 0, otherwise. 

The limiting behavior of the spectrum of T n (A s ) as n — > oo is characterized 
in terms of the solutions of the algebraic equation A s (z) — x, see [2J. For every 
x G C there exist exactly m — 1 solutions of the equation A s (z) = a; (assume that 
6( m ~ 2 )(s) ^ 0) which we denote by zj(x, s), j = 1,.. . , m — 1. We label the solutions 
by their absolute value so that 

(1.15) \zi{x,s)\ > \z 2 (x,s)\ >■■■> \z m -i{x,s)\ > 0. 
We put 

(1.16) r x (s) = {xeC \ \zi(x, s)\ = \z 2 (x, s)\}, 

which is a finite union of analytic arcs. 

We will use the following classical theorem on the behavior of the eigenvalues of 
T n (A s ) as n — > oo. 

Theorem 1.1. As n — > oo the eigenvalues of T n (A s ) accumulate on the contour 

(fTTB) . 

The sequence of normalized counting measures of the eigenvalues ofT n (A s ) con- 
verges weakly as n — > oo to the Borel probability measure [i\ on Yi(s), given by 

1 ( z[_(x, s) z' 1+ (x, s) 



(1.17) dri(x) = — - -^-^ dx, 

2m \zi-(x,s) zi + (x,s) J 

where ' denotes the derivative with respect to x. In (|1.17|) . we have that dx is 
the complex line element on Ti(s) and Zi±(x,s) is the limiting value of zi(x,s) 
as x — * x from the ± side on each of the arcs in Ti(s) (the complex line element 
induces an orientation on ^i(s) and the + side (— side) is on the left (right) as 
one traverses Y\(s) according to the orientaton). 

Proof. The fact that the eigenvalues accumulate on Ti(s) was shown by Schmidt 
and Spitzer |16j . The result about the limit of the normalized eigenvalue counting 
measures is due to Hirschman [10], see also [2j Chapter 11]. The precise form (|1.17|) 
for the limiting measure was given in [6]. □ 

For later use we also note that by |3 Proposition 4.2], 

(1.18) fMM = fi^4 fo r, e C\r l(s ), 
J x-y zi[x,s) 

which also characterizes the measure /xf . 

Our first result states that, under certain conditions, the polynomials Pk iTl sat- 
isfying the recurrence (JTTT2J) have a limiting zero distribution, which is an average 
of the measures (jl.171) . The average is with respect to the parameter s. 

Theorem 11.21 is an extension of Theorem 3.1 of [3] to polynomials satisfying an 
m-term recurrence relation instead of a specific four-term recurrence relation as 
in [3]. The analogous result for orthogonal polynomials satisfying a three-term 
recurrence is from [14]. See also [13] and [17] for related results. 
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Theorem 1.2. Let for each n e N, m — 1 sequences {&jjf'^}jfcLo> J = 0, . . . , ttz — 2, 
o/ reaZ coefficients be given and assume that there exist continuous functions b^' : 
[0, oo) — * R, j = 0, . . . , m — 2, suc/i i/ioi /or eac/i s > 0, 

(1.19) Um bf n = 6 (i) (s), j = 0, . . . , m - 2. 

£e£ -Pfc,n &e the monic polynomials generated by the recurrence (|1.12|) and suppose 
that 

(a) £/ie polynomials Pk,n have real and simple zeros x\ ,n < ■ ■ ■ < x k ' n satisfying 
for each k and n the interlacing property 

Xj < Xj' <x j+1 ' , for j = 1, . . . ,k, 

(b) Ti(s) C R /or every s > 0, where Ti(s) is given by ()1.16p . 

TTien i/ie normalized zero counting measures v{Pk.n) — j Ej=i /lave a weafc 

Zirmi as fc, n — » oo wzi/i fc/n — ► £ > given by 

(1.20) lim j/(P fe ,„) = - / /i^ds 

fc/n->£ t, Jo 

where /xf is the measure (I1.17[) . 

The proof of Theorem 11.21 is given in Section [2] 

1.4. Multiple orthogonal polynomials associated with modified Bessel 
functions. We want to apply Theorem 1 1.21 to the multiple orthogonal polynomials 
Bk associated with the modified Bessel function. 

1.4.1. Interlacing. The assumption (a) of Theorem 11.21 will be satisfied since we 
have the following. 

Proposition 1.3. Let a > 0, a > —1, < t < T. Then the polynomials Bk 
generated by the recurrence (|1.8p with recurrence coefficients (|1.9p have real and 
simple zeros in (0, oo) with the interlacing property. 

The proof of Proposition 11.31 is given in Section [3] 

1.4.2. First reseating. In the first rescaling we let t and T depend on n, while a 
and a remains fixed. We replace 

t^t/(2n), T^l/(2n) 

and we obtain the recurrence coefficients as in (jl.lljl . The scaling limits of the 
recurrence coefficients indeed exist: 

lim b k .n = b(s) = a (l - tf + 2st i l - *)> 

k/n — >s 

(121) I™ c k .n =c(s) = 2ast{l-tf + s 2 t 2 {l-t) 2 1 

\ ' k/n^s 

lim d k<n = d(s) = as 2 t 2 (l - t) 4 . 

k/n — >s 

Then as in (|1.13p we have the associated family of symbols 
(1.22) A s {z) = z + b(s) + c(s)z" 1 + d{s)z- 2 
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Figure 2. Graph of the curves /3(s), j(s) and rj(s) for the first 
and second rescaling. 



and the solutions z\(x, s), z-z(x, s) and Zs(x, s) of the algebraic equation A s (z) = x. 
We define T\(s) as in (|1.16p and similarly 

(1.23) r 2 (s) = {x G C | \z 2 {x, s)\ = \z 3 {x, s)\}. 

The symbol (|1.22p with the functions b(s), c(s) and d(s) from (|1.21[) allows for 
a factorization 

(1.24) {z + a{l-tf){z + st{l-t)f_ 

z 1 

From (|1.24[) we see that A s has three negative zeros, namely a double zero at 
—st(l — t) and a simple zero at — a(l — i) 2 . For the special value 

* a(l-t) 

s = s ' = 

i 

the three zeros of the symbol coincide. 

These facts are used to prove the following. 

Proposition 1.4. For each s > 0, we have that Ti(s) C [0, oo) and ^(s) C 
(—oo,0]. More precisely, there exist r)(s) < < /3(s) < 7(5) so i/iai 

(1.25) r 1 ( s ) = [/3( s ), 7 ( s )], r 2 ( S ) = (-oo,7j( S )]. 

In addition we have 

(a) s 1— > 7(5) is strictly increasing for s > 0, lim s ^ + j(s) — a(l — t) 2 , and 
lim^oo 7(s) = +oo ; 

(b) s 1 — ► /?(s) is positive and strictly decreasing for < s < s* = "(i - *) and 
/3(s) = for s > s* . Furthermore, lim s ^ + P( s ) — a (l — t) 2 , 

(c) ry(s) = for < s < s* and s > r?(s) zs negative, strictly decreasing for 
s > s* and lim^oo r;(s) = —00. 

From the proposition it follows that the sets Ti(s) and ^(s) are intervals, and 
that Ti(s) is increasing as s increases, while ^(s) decreases. See the left part of 
Figure [2] for an illustrative plot of 'j(s), (3(s) and T)(s), as a function of s. 
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As a result of Propositions 11.31 and 11.41 we see that the assumptions (a) and (b) 
of Theorem [T72] are satisfied, and so for each £ > 0, the weak limit of the normalized 
zero counting measures of the polynomials as k,n — > oo, k/n — > £, exists and 
is given by (|1.20[) . 

1.4.3. Second rescaling. In the second rescaling we let the parameter a increase as 
n increases. For this we change the variables t, T and a by 

t 1 ► t J (2ti) , T \- >l/(2n), ai-tpn 

with p > 0. 

The recurrence coefficients (|1.9j) now have scaling limits that also depend on p. 
Indeed, 

b(s) = a(l - t) 2 + 2st(l - t) + t(l - t)p, 
(1.26) c(s) = 2as<(l - t) 3 + s 2 t 2 (l - t) 2 + st 2 (l - t) 2 p, 

d(s) = as 2 t 2 {l-tf. 
Now the symbol (11.24)) depends on p and we can again explicitly factorize 

(127) AJz) - {Z + St{1 - t)){z2 + (1 ~ t)(a(1 - t] + ( ' S + aSt{l - t)3) 

* z 2 

There are again three negative zeros of the symbol, but now all three zeros are 
simple. 

We again have Ti (s) and T2 (s) as in (|1.16[) and (|1.23p , and we prove the following. 

Proposition 1.5. Let p > 0. Then for each s > 0, we have that Ti(s) C (0,oo) 

and ^2(s) C (—00, —^7—]. More precisely, there exist r)(s) < < (3(s) < j(s) so 
that 

(1.28) T^s) = [(3(s), 7 («)], T 2 (s) = (-00, r,(s)}. 

In addition we have 

(a) s 1 — ► 7(s) is strictly increasing for s > 0, lini s ^o + 7(s) = (1 — t)(a(l — t)+pt) 
and lims^oo "/(s) — 00, 

(b) s 1 — ► /?(s) is positive and strictly decreasing for s > u;ii/i lim s ^o+ = 
(1 - f)(a(l - i) +pt) and lim^oo [3{s) = 0, 

(c) s 1 — ► 77(5) is negative and strictly decreasing for s > wii/i lim s ^o+ 7(s) = 

and lim^oo ry(s) = -00. 

See the right part of Figure [2] for a plot of the functions (3(s), 7(5) and n(s) in 
the case p = 5. 

Also in the second scaling the assumptions of Theorem 11.21 are satisfied, and 
again it follows that for each £ > 0, the weak limit of the normalized zero counting 
measures of the polynomials Bk, n as k, n — ► 00, k/n — > £, exists and is given by 
(QUI) . 

The proof of Propositions 11.41 and 11.51 are given in Section 01 

1.5. Equilibrium problem. In both scalings we find for each £ > a probability 
measure of the form 

1 r 6 

(1.29) "1 = 7 Vids 

£ Jo 
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as the weak limit of the normalized zero counting measures. The main result of the 
paper is that this measure can also be obtained as the first component of a vector 
of measures (vf, v 2 ) that satisfies a vector equilibrium problem. For the case p = 
and £ = 1, it reduces to the vector equilibrium problem (|l.l|) - (|1.4|l stated in the 
introduction. 

We will use a recent result of Duits and Kuijlaars [B] which in the present context 
says that the measure /zf that gives the limiting eigenvalue distribution of the 
Toeplitz matrices with symbol (jl.24l) is part of a vector (/xf , fi 2 ) that is characterized 
by a vector equilibrium problem. 

The second measure \i\ is supported on T 2 (s) (see (|1.23[l ) and is given by 

(1.30) dfi 2 (x) = — ±- dx 

2m \z 2 -(x, s) z 2+ {x,s)J 

for x G r 2 (.s). Then fj, 2 is indeed a positive measure on T 2 (s) with total mass 1/2. 
The result of [Sj in this special case is the following. 

Theorem 1.6. For each s > we have that the vector (A*f,/x|) is the unique 
minimizer for the energy functional 

(1.31) //log— ^ — | dm(x)dfjLx(y) + //log I 1 , dfi 2 (x)dfi 2 (y) 

JJ \%-y\ JJ \x-y\ 

log- -dm(x)dfj, 2 (y) 

F - y\ 

among all vectors (^i,^l 2 ) satisfying supp(/Zj) C Tj(s) for j — 1,2, and 

dm = 1, J dp, 2 = ^. 
The measures and fi 2 satisfy for some constant I s , 



(1.32) 2J \og\x - y\dnl(y) - J log\x - y\dfj, s 2 (y) = t , xGT^s), 

(1.33) 2 J \og\x - y\d^ 2 (y) - J \og\x - y\dfil(y) = 0, x £ T 2 (s). 

The conditions (|1.32|) - (|1.33|l are the Euler-Lagrange variational conditions for 
the vector equilibrium problem. 

Recall that we have (|1.29l) and similarly we put for £ > 0, 

c 1 ft 

(1.34) 4 = 7 ^ ds - 

s Jo 

Then v\ is a measure on U s <f ^(s) = (— oo, — p 2 i 2 /4a] with total mass 1/2. We 
obtain the vector equilibrium problem for {v\, v 2 ) by integrating the vector equilib- 
rium problem for (/if, with respect to s, in particular the variational conditions 
(ll.32p - (|1.33p . A complication is that the intervals T\(s) and T 2 (s) are varying with 
s. The fact that Ti(s) is increasing with s induces, after integration, an external 
field on vf. The fact that r 2 (s) is decreasing as s increases leads to the upper 



constraint on v 2 



2 ■ 

Theorem 1.7. Define 



(1.35) V(x) = / log 



zi{x,s) 



z 2 (x, s) 



ds 
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a=l,t= 0.2, p = Q, £_ = 1 
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X 

FIGURE 3. Graphs of the densities of a (dashed) and v 2 on the 
negative half line and the density of v\ on the positive half line, 
for the case p = and t < t* . 



and 

poo 

(1.36) a = I n s 2 ds. 

Jo 

Then for every £ > 0, the vector of measures (^f,^!) * s the unique minimizer 
for the energy functional 

(1.37) / / \og—^—d Vl {x)d Vl {y) +[[ \ og ——dv 2 {x)dv 2 (y) 
J J \x-y\ J J \x-y\ 

1 If 

log- -dv 1 (x)dv 2 {y) + - / V(x)dv 1 (x), 

\x-y\ U 

over all vectors of measures (^1,^2) such that supp(^i) C [0,oo), J dv\ = 1, 
supp^) C (— oo,0] ; J dv 2 — 1/2, and 

1 

v 2 < -a. 

The measures v\ and v\ are characterized by the following variational conditions 
(1.38) 

2/loglx-^fa,)- [log\x -y\d4(y) + jV(x)H; f OT ^™^)' 
J J ? [< t, for x € [0,oo), 

for some £, and 
(1.39) 

1 /"i 1 1,7 if \ /"i 1 1, ^ J=°. for a; G supp (o--^|) ; 

2 / log \x-y\dvl(y - / log \x-y\dv\(y) I * 

J J l>0, for x e C \ supp(cr - 

The proof of Theorem 11.71 is given in Section [5] 
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o= 1, t = 0.9, p = Q, f = 1 












" _jA_^y 





-1.6 1 -0.5 0.5 

a; 



Figure 4. Graphs of the densities of <r (dashed) and j/f on the 
negative half line and the density of v\ on the positive half line, 
for the case p — and t > t* . 



1.6. Evaluation of V and er. In a final result we are able to evaluate the integrals 
(fj~35f and (fl~36)l that define V and a. 

Theorem 1.8. For every p > 0, we /icroe 



(1.40) = - ^ pH2 f + - plog (Vp 2 * 2 + 4a* - pt) 



£<l-2+plog(2a(l-t)), 



and a is the measure on (— oo, 0] wi/i density 
da(x) 



(1.41) 



\/4ojxf— p 2 t 2 
27rtH ' 



fix 








Note that for p = 0, the external field (|1.40p is 

<" 2 » ^' = j(i^)-^ + ^ = j(rh)(^-^ 1 - i »» 1 

and the constraint p.4ip is the measure with density 

|xr 1/2 , a? G (-oo.O]. 



dcr(a;) _ v / a|„i-i/2 



dx irt 

Remark 1.9. Figure [3] shows the graph of the densities of a, v\ and v\ in the case 
p = and i below the critical time i*. This corresponds to the case where the 
non-intersecting paths did not come to the hard edge. The constraint a is not 
active and v\ is supported on a interval which is at a positive distance from zero. 
Figure [4] illustrates the case p = for t > t* . Here the constraint a is active in an 
interval. 

Figure [5] shows the densities of a, vf and v\ for p > 0. For all values of t, the 
constraint a is active in some interval and v\ is supported on an interval which is 
at a positive distance from zero. 
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1 

a = l,t = 0.9, p = 1, £ = 1 




i i 












V 


i i i 






\ 

i i 1 i 



-3 -2-10 1 2 3 

FIGURE 5. Graphs of the densities of a (dashed) and v\ on the 
negative real line and the density of v\ on the positive real line for 
p > 0. 



Remark 1.10. It was pointed out that if we let n — > oo, the paths fill out a region 
in the te-plane. This can be observed in the left figure of Figure Q] for the case 
p = and in the right figure of Figure Q] for the case p > 0. The region filled by 
the paths, for a fixed time t 6 [0,1], is exactly the interval [(3(1), 7(1)]- We can 
obtain (3(1) and 7(1) by computing the zeros of the discriminant of the polynomial 
z 1 A\(z) — z 2 x with respect to z, which is the following algebraic equation of degree 
three in x: 

Aax 3 - [8a 2 (l - tf + Aat(l - t)(2p + 5) - t 2 (p + lf]x 2 
+ (l-t)[4a 3 (l-t) 3 +^a 2 t(l-t) 2 (2p-3) + 2at 2 (l-t)(p 2 +p+6)-2t 3 (p+2)(p+l) 2 ]x 
+ p 2 t 2 (l - t) 2 [a 2 (l - t) 2 + 2at(l -t)(p-l)+ t 2 (p + l) 2 ] = 0. 

If p = then the algebraic equation reduces to 

(1.43) x(4ax 2 - ir(8a 2 (l - t) 2 + 20at(l - t) - t 2 ) + 4(1 - t)(a(l - t) - t) 3 ) = 0. 
We can compute explicit expressions for the solutions in this case 
xi(t) = 0, 

x 2 (t) = ^~(8a 2 (l - t) 2 - t(t - 20a(l - t)) - y/t(t + 8a(l - t)) 3 ) , 

x 3 (t) = — (8a 2 (l - t) 2 - t(t + 20a(l - t)) + yjt(t + 8a(l - t)) 3 ) . 
8a 

Remark 1.11. If we let a — * and s = 1, then the symbol (|1.27|) becomes 
(l + zt(l~t)) 2 pt(l-t)(l + zt(l-t)) 
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1.2 




Figure 6. Numerical simulation of 50 non-intersecting paths for 
a = and p = (left), p > (right). 



The solution z 3 (x, 1) tends to zero as a — > 0. On the other hand, the solutions 
zi(x, 1) and z 2 (x, 1) have limits (with appropriate choice of ±-sign) 

- ~ *)(? + 2) + \ + \ ((* - Pi(t))(* - P2(<))) 1/2 , 

*a(*) = -^(1 - *)(? + 2) + \ ~ \ ((* ~ Pi (*))(* - P 2 (<))) 1/2 , 
as a — > 0, where 



pi(t) =i(l-t)(p + 2-2Vj»+l), p 2 («) = i(l - t)(p + 2 + 2Vp + 1). 
One can show that v\ is the Marchenko-Pastur distribution, see e.g. [S], with density 



rf^i(x) v / (p 2 (Q-x)(a;-pi(f)) 

-*T = M(l-t)x ' ase ^W'^*)]- 

Figure O shows simulations of 50 non-intersecting squared Bessel paths for the case 
a = and the boundaries yOi(t) and p 2 (i) of the region filled by the paths as n — > oo. 

The rest of the paper is devoted to the proofs of the theorems and propositions 
stated above. Theorem 11.21 is proved in Section [51 Proposition 11.31 in Section [31 
Propositions 11.41 and 11.51 in Section [4] Theorem 11.71 in Section and finally Theo- 
rem 11.81 is proved in Section [51 



2. Proof of Theorem 11.21 

Before proving Theorem 11.21 we will state a result concerning ratio asymptotics 
of polynomials satisfying a recurrence relation (|1.12p with varying recurrence coef- 
ficients. 

The following lemma will be used in the proofs of Theorem 11.21 and Lemma l2~2l 

Lemma 2.1. Suppose that the zeros of the monic polynomials Pk and Pk+i, with 
degrees k and k + 1, respectively, are simple and real, lie in an interval [-R, R] for 
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some R > 0, and are interlacing. Then we have 



(2.1) 
(2.2) 

(2.3) 



Pk(z) 



Pk+i(z) 
Pk(z) 



< 



1 



dist(z, [-R,R]) 



forze C\[-R,R], 



Pk{z) 
Pk+i(z) 



> ^ for \z\ > R, 



< 



dist(z, [-R,R]) 2 



forzeC\[-R, R]. 



Proof. Inequalities (|2.ip and (|2.2p can be found in []3J Lemma 2.2]. The proof of 
is similar to that of (|2.1[) . □ 



We consider the doubly indexed sequences of polynomials {Pk,n\ generated by 
the m-term recurrence relation (|1.12|) and we assume that the recurrence coefficients 
have scaling limits 



lim b 



k/r. 



U) 

k,n 



Lemma 2.2. Under the assumptions of Theorem ] 1.2\ we have that for each s > 0, 
there exists R > so that all zeros of Pk, n belong to [-R, R) whenever k < (s + l)n. 
Moreover, 



(2.4) 



lim fc+i-,"^) = Zi(x, s), 



k/n^s P k ,n(x) 

uniformly on compact subsets ofC \ [—R,R]. 

Proof. Fix s > 0. The convergence (|1.19[) imply that the recurrence coefficients 
are uniformly bounded if k/n is restricted to compact subsets of [0,oo). So, the 
number R defined by 

(2.5) R := sup{l + |6gj + + ■ • • + \b { ^ 2) \ \ k < (s + 1)„} 

is finite. 

By the recurrence (|1.12l) we have Pk, n {x) — det(zlf- — Mk, n ) where M/, y n is the 
matrix 



/ 6 (0) 

/ "0.71 



(2.6) 



,(m-2) 
J m-2,n 





1 

l(0) 



V 



o 



,(m-2) 
-'fe-l.n 



1 

do) 



b W / 



The zeros of Pfc jrl are equal to the eigenvalues of Mk, n > an d their absolute values are 
bounded by the maximum absolute row sum of Mk, n - Therefore, by the definition 
3)) , if k < (s + l)n then the zeros of Pk,n lie in the interval [-R, R]. 
We consider the family of functions 

ffc+l,?l 



(2.7) 



n = 



| k,n e N, k + 1 < (s + l)n 
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Because of the assumption (a) in Theorem ll.2l we can apply Lemma I27T1 to Pfe jn and 
Pk+i,n- It follows from (12.11) . (|2.2p and (12. 7|) that the family Ti is a normal family 
(in the sense of Montel) on C \ [— R, R}. 

Using induction on I, we will show the following. 

Claim: For each I > 0, the following holds. If {ki}i, {rii}i are sequences of 
non-negative integers with ki,rii — » oo, kijrii — > s as i — > oo, so that 

-Pfci+i,n 4 (a;) 



(2.8) 



lim 



exists for |x| > R, then 
(2.9) f(x) = zi(x,s)+0(x- 1 ) 

as x — > oo. 

We have zi(x, s) = x + 0(1), as a? — > oo, and so it is clear that (|2.9[) holds for 
1 = 0. 

Now assume that the claim holds for I > 0. Let {fc^} and {ni} be as in the claim. 
Since kijrii — > s as i — > oo, we may assume that fcj < (s + l)rii for every i. For 
j = 0, . . . , m — 2, we then have that 

Pki + l-j,m 

belongs to the family 7i. Since it is a normal family, we may assume, by passing to 
a subsequence if necessary, that 



(2.10) 



/ (i) (x) = lim 



-oo P k ._ jn .(x) 



exists for x E C \ [-R, R) and j = 0, . . . , to — 2. 

Now we divide the recurrence (|1.12[) by Pfe in , and replace k and n by ki and rij, 
to obtain for each j, 



-Pfci+i,nj(a:) , (0) 



E 6 



0) ffc,-j>,(x) 



We let i — > oo, where we note that by (|1.19p 



1,(3) 



6«(s) 



and by (j2~Tu| 



/«(x)---/W(x) 



for x S C \ [-i?, i?] and j = 1, . . . , to - 2. Thus 

m-2 

x = /(x) + (s) + bU) ( s ) ( x ) ' ' ' f U) ( x ) 

3 = 1 

for all x € C\ 

Now by the induction hypothesis we have for j = 0, . . . , to — 2, 



= Zl (x, s) + 0(x-') = Zl (x, s) (1 + ©(x^ 1 )) 



as x — ► oo. 
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and so 

m-2 

f(x) = x- &(°) (s) - £ 00 (?) ■ ■ ■ f U) (*) 

3 = 1 
m-2 

^x-b^is) - b^{s)z 1 (x,s)- : > (l + Oix- 1 - 1 )) 

m-2 

= x-b^(s)- b^\s)z 1 (x,s)- j +0(x~ 1 - 2 ), 
i=i 

Since zi(x, s) is a solution of = x, we have 

m-2 

a;-6 (0) (s)- b U) {s)z 1 {x,s)- 1 =z 1 (x,s) 

and so we obtain 



f(x)=z 1 (x,s)+0(x- 1 - 2 ) 



as x — * oo, which proves (|2.9[) for i + 2. 

The claim now proved, we finally show how the lemma follows from the claim. 
First note that / is analytic in C\ [— R, R] and x t— > zi(ir, s) is defined, analytic and 
non-zero in C \ Thus, if (|2.9[) holds for every I, then clearly f(x) = Z\(x, s) 

for a; in a neighborhood of oo in the x-plane, and by analyticity the equality extends 
to C\ ([— R, R]UTi(s)). Next, recall the assumption (b) in Thcorcm ll.21 which says 
that Ti(s) C R. Then it easily follows by its definition that Z\{x, s) cannot possibly 
have an analytic continuation from C\Ti(s) to a larger set. Thus T\(s) C [—R,R] 
and it follows that 

(2.11) Urn P fe«+b"iW = Zl ( XjS ), xeC\[-R,R], 

woo Pk^nAX) 

for all sequences {fci}, as in the statement of (|2.7p . 

Then by a standard normal families argument (recall that the family TL is a 
normal family) (|2.4p follows. 

□ 



Proof of Theorem 11.21 Let £ > 0, and take < s < 1. From Lemma 12.21 there 
exists R > such that all zeros of Pk, n belong to [—R,R]. For k £ N, we use |_s&J 
to denote the greatest integer less than or equal to sk. We observe that 



kPk,n(x) k \Pj +ltn (x) Pj,n(x) 

(2.12) = / ffi^^ - ff'Hl W 



By taking the logarithmic derivative of (|2.4|) and using (j 1 . 18[) we obtain 

/ 213 s Um f P [sk}+l,n( X ) P jsk\,n( X ) \ = Z[(X,S£) _ f dfxf (y) 
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uniformly on compact subsets of C \ [-R, R]. From (12. 2\ and (|2.3| in Lemma |2~T1 
we obtain 



(2.14) 



Ysk\, 



LsfeJ + l.n^) P[sk\,n(x) 



Pi 



< 



[sk]+l,n 

2\z\ 



(x) 



P[sk] +l,n(x) 



Pi 



dist(z, [-R,R]) 2 

for |2| > R. Therefore we may apply Lebesgue's dominated convergence theorem 
and it follows by (|2TT2"j) and (pTTB|> 



lim 



1 ^nW 

Therefore we have that 

lim 



ds 



^Mi(y) 



ds. 



ds. 



k/n->£ k Pk,n{x) 

for all x E C \ [— R, R]. This gives, by a standard argument, see 
polynomials Pk. n have 



Mi(A)ds, 



as limiting zero distribution as k/n — > £. 

3. Proof of Proposition LL3I 



that the 



□ 



In [4] it is proved that the weights (wi , w 2 ) from (|1.5p are such that there exist 
discrete measures <j\ and a 2 on (— oo,0] such that 



w 2 (x) 
wi(x) 



d<j\{i) wi(x) 



da 2 (t) 



x>0. 



, x — t ' 102(2;) 

From this it is shown in [4j Theorem 4] that the weights (w±, w 2 ) form an AT-system 
on [0, 00), which means that for any m, n 2 € N, the set of real valued functions 



J n 



{wi(x),XWl(x), . . . ,X ni Wl(x),W 2 (x),XW 2 (x), 



,x n2 w 2 (x)}, 

U1+U2+2 

I 

0), has at 



is a Chebyshev system on (0, 00) i.e., every linear combination Yuk=d ~ a k<Pk, 
with ip k £ Fnunv Vk ^ <Pk' if ^ k', and (ai, . . . , a n+m+2 ) ^ (0, . . . , 
most m + n 2 + 1 zeros in (0, 00). 

A consequence of this result is that for every n\ , n 2 £ N there is a unique monic 
polynomial B niiri2 of degree ni + n 2 satisfying the orthogonality conditions 



(3.1) 



B ni ,n 2 (x)x k Wj(x)dx = 0, 



for k = 0,1,. 



3 = 1,2. 



The polynomial B ni ^ n2 has exactly rti + ^2 simple zeros in (0,oo), see e.g. |19) . 
Comparing (|3.ip with (|1.5|) we see that B n = B niTl2 with n\ = n 2 = n/2 if n is 
even, and ni = (n + l)/2, n 2 = (n — l)/2 if n is odd. 

The interlacing of zeros now follows from Lemma 2.3 and Remark 2.1 of pQ. 
Indeed there it is shown that for any m, n 2 the zeros of B nitTl2 are interlaced with 
those of B ni+ i >rL2 and with those of B ni _ rL2+ i. 

For convenience of the reader, we give the proof of interlacing following p] (see 
also [H Lemma 5 and Corollary 1] for a different approach). We show that the zeros 
of B nit n 2 interlace with those of B Ul +i,n 2 , the proof for B ni ,n 2 +\ being similar. 
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/ 7 ' 

/ H 

^s. V ? J 


riW 




fr\ y'i r$\. 


A3 

1 


2/3 2 

r 2 (s) 



FIGURE 7. The graph of A s {z) (fL2"Tj) . for p > 0. 



Let us consider a polynomial P — aB ni7l2 + bB„ 1+ i tn2 with (a, 6) 7^ (0,0). Let 
us assume that P has a double real zero at ( £ K. Then 

(3.2) P(x) = aB nun2 (x) + bB ni+1 , n2 (x) = {x-() 2 R{x), degR < n x + n 2 - 1. 
From the orthogonality conditions (|3.1|) it follows that R satisfies 



(3.3) 



R(x)x k (x - () 2 Wj(x)dx = I (aB nitri2 (x) + bB ni+1 . n2 {x)) Wj(x)dx 



for k = 0,...,nj - I, j = 1,2. 



Thus i? has multiple orthogonality conditions with respect to the weights (wi,w 2 ) 
where 



wi(x) = (x - C) 2 wi(x), 



and 



w 2 (x) = (x- () 2 w 2 (x). 



By the same reasoning as in [3] it follows that (wi,w 2 ) is also an AT system on 
[0,oo). It follows that any monic polynomial R that satisfies (|3.3[) should have 
degree > n\ + n 2 , and this is a contradiction. 

Therefore the polynomial P has only simple zeros in R. This means that for 
every fixed x € R, the linear system 



/ a 

B' ni , n2 (x) B' ni+hn2 (x)J \b 



has only the trivial solution a = b = 0. Then the determinant is non-zero 

B ni . n2 (x)B' ni+hn2 (x) - B ni+hn2 (x)B' nin2 (x) 7^ 0, 

which readily implies that B nitTl2 has opposite signs at consecutive zeros of i? ni +i.„ 2 . 
It follows that between two consecutive zeros of B ni+ \^ n2 there is at least one zero 
of B niv n 2 , and this proves the interlacing of zeros. □ 
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4. Proof of Propositions IL4I and LL5 



Propositions 11.41 and 11.51 will be proved simultaneously. First of all we observe 
that the symbol allows for a factorization 

(4.1) A s{z) = {Z ' ri){Z \ r2){Z ' n \ with n, r 2 ,r 3 <0, 

z z 

see (|1.24j) for the case p = 0, and (j 1 . 2T|) . We order the roots so that 

ri < r 2 < r 3 < 0. 

If p > 0, then it is elementary to check from (|1.27p that all zeros are distinct. 
Figure [7] shows the plot of A s (z) for the case p > 0. If p = 0, then it follows from 
(fj~24f that 

ri = r 2 < r 3 < 0, if st < a(l - t), 
n = 7'2 = r 3 < 0, if st = a(l — t), 
n < r 2 = r 3 < 0, if st > a(l - t). 

These three cases are illustrated in Figures 13 [9] and [HH respectively. 
The derivative of A s (z), 

A' s (z) = 1 - c{s)z- 2 - 2d(s)z~ 3 , 

has three roots in the complex plane. From Figures [3 EJ El and fTOl we see that all 
zeros are of A' s are real. We denote the zeros of A'(z) by yi, y% and y 3 so that 

Vi < V2 < < y 3 , 

as indicated in the figures. The equality y\ = yi only holds if p = and s = s* as 
in Figure [9l To indicate the dependence on s we also write yi(s), j/2(s) and y 3 (s). 

Before proving Propositions 11.41 and 11.51 we first need two lemmas. The proofs 
of these lemmas only use the fact that the zeros rj of the symbol A s are strictly 
negative. 

Lemma 4.1. Assume that Zi,z% G C are such that Z\ ^ zi, \z\\ — \zi\ and 
A s {z\) — A s (z2) — x. Then z± — zi and 

Proof. The complex numbers z\ and zi lie in a circle of radius p — \zi\ — \z2\ 
centered at the origin in the complex plane. Then from (|4.ip we have 

dist(2,n)dist(z,r2)dist(z,r3) . . 

\A s {z)\ = , a\z\=p. 

P 

Since all rj < 0, it follows that 

[— 7T, 7t] — > R : 8^ \A s {pe w )\ 

is an even function which is strictly decreasing as increases from to n. 
Thus equality 

\A s (pe^)\ = \A s (pe ie *)\, 
with #1.2 S [— 7r, 7r] and £q 7^ 2 can only occur if 2 = Then it follows from 

the assumptions of the lemma, that z\ = Z2, and 



x = A s {z 1 ) = A s (zi) = A s {z 2 ) = x 
so that x GR. □ 
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A s (z) 


7 v_ 

ri(«) 
7? 








r 2 ( S ) 



Figure 8. Graph of A s (z) (|1.24|) . z e 1, in the case p = 0, s < s* . 

Lemma 4.2. For each s > we have Ti(s) U ^(s) C R. Moreover, if p > or i/ 
p = and s 7^ s* = a(l — £)/£, i/ien Ti(s) n ^(s) = 0. On the other hand, if p = 
and s = s* , then Fi(s) fl T2(s) = {0}. 

Proof. If x G ri(s)ur2(s) and — x has a double root, then there exists z\ G C 

such that A s (zi) = a; and v4' s (zi) = 0. Since all zeros of A'(z) are real, it follows 
that zi G K. 

If a; G Ti(s) U T2(s) and A s (z) — x does not have a double root, then there exist 
z\,z 2 G C such that Z\ ^ z 2 , \zi\ = \z%\ and x — A s {z\) = A s (z2)- Lemma |4~T1 savs 
that x G M. This proves that T^s) U r 2 (s) C M. 

If x G ri(s) n r 2 (s), and s ^ s* if p — 0, then there exist three solutions of 
J 4 s (z) = x. One negative solution Z\ < and two complex conjugated solutions 
2:2 and Z2 (see Figures [3 [lOj) . Moreover z\ ^ Z2 and z\ ^ z 2 - Now |zi| = |z 2 | 
(a; G ri(s)) and A s {z\) = A s {z2). Therefore we can apply Lemma 14.11 and we 
obtain z\ = z 2 which is a contradiction. Then ri(s) n T2(s) = 0. 

The fact that Ti(a*) n r 2 (s*) = {0}, if p = 0, follows directly from Figure El by 
observing that z = —a(l — t) 2 is a triple root of A s (z). □ 

Proof of Propositions 11.41 and 11.51 If p > then there are three local extrema of 
A s (z), namely f3(s),j(s),rj(s), such that 

n(s) < < 0(s) < 7(5). 

If x € (rj(s), f3(s)) U (7(5), 00), then there exist three different real solutions of 
A s (z) = x. It is easily seen that these solutions also differ in absolute value (Fig- 
ure[71 and Lemma l4.1l if x G (7(5), 00)). On the other hand, there is one real and 
two complex conjugated solutions whenever x G (—00, rj(s)] U [/3(s), 7(s)] . Therefore 

r x (s) U r 2 (s) c (-00, ri(s)] U \p(8)M*)]- 

We have that A s (z) — j(s) has a double root y 3 = y 3 (s) > and one negative 
root whose absolute value is less than 2/3 (s). Thus 7(s) G ri(s). Also /3(s) G ri(s). 
Now the fact that Ti(s) is connected (see [H],[21 Theorem 11.19]) says that Ti(s) — 
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FIGURE 9. Graph of A s (z) (|1.24|) . z e 1, in the case p = 0, s = s*. 



[/3(s),7(s)]. On the other hand A s (z) = T](s) has a double root at 2/2(5) < and a 
negative root whose absolute value is larger than 1 2/2 C s ) I - Therefore 772(5) € ^(s) 
and T 2 (s) = (-00, rj(s)]. 

The case p = follows analogously by studying Figures El El and [TO] whenever 
st < a(l — t), st = a(l — t) or st > a(l — t), respectively. We deduce from Figure [5] 
that rj(s) = for st < a(l — t). On the other hand, Figure [TO] shows that /3(s) = 
if st > a(l - t). 

We shall prove that 7(s) is an increasing function. The monotonicity of (3(s) and 
rj{s) is proved with similar considerations. Observe that the function 

_ (££ + a(l - tf){z + t(l - t)f pt(l - t){z + t(l - t)) 

-D{Z, Sj — Si s (SZ) — 7. 1 , 

z z z 
as a function of z, has a local minimum at j/3(s)/s and 

7(a)=^.(l/s(*)) =B(y 3 (a)/a,s). 

If we take the partial derivative of -B(z, s) with respect to s we obtain 

dB(z,s) _ (z + t(l-t)) 2 
ds z 

which is positive for all z > 0. Then the fact that z t— > B{z, s) has a minimum at 
ys,{s)/s implies that 

, . = dB(y 3 (s)/s, s) (y' 3 {s)s - y 3 (s)) + dB{y 3 {s)/s, s) = dB(y 3 {s)/s,s) > Q 
dz s 2 ds ds 

and therefore j(s) is increasing for all s > 0. 
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FIGURE 10. Plots of A g (z) (fOij) . z e K, in the case p = 0, 
s > <z(l- 

If p > 0, it is a straightforward computation that yi(s), 2/2(5) and 1/3(5) have the 
following behavior as s — > 00 

yi (s) = st(l-t) + 0(l), 

y 2 ( S ) = -2a(l-t) 2 + 0( S - 1 ), 

y 3 (s) = -st(l - t) - |t(l - i)p + ©(s- 1 ). 

Then we obtain 

7(a) =^(lfc(«))=45f(l-t) + 0(1), 

/?(*) = MM*)) = fp 2 ( l - + o{s~ 2 ), 
v ( s ) = Myi(s)) = -^ + o(s), 

as s — > 00. This proves the limits 

lim /3(s) — > 0, lim 7(s) = 00, lim r}(s) = —00, 

s — >oo S — >00 S — >00 

and completes the proof of the propositions. □ 

5. Proof of Theorem 11.71 

From the definitions of v\ and v\ in (|1.29[) and (|1.34p , it is immediate that 
supp(z^) C [0,oo), J du\ — 1, supp(z/f) C (—oo,0] and J dv^ = 1/2. In fact, from 
(|1.29p and the fact that the sets T\(s) — supp(/xf) are increasing as s increases (see 
Propositions 11.41 and 1 1 . 5 j) it follows that 

(5.1) su PP (4) = |Jr 1 ( s ) = r 1 (o. 
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From (|1 .34|) and the definition of a in (|1.36j) . it is then also clear that 



A - 



in i r°° i 

- J fi s 2 ds < - J fi s 2 ds = -a, 



and 
(5.2) 

so that 
(5.3) 



a - £v\ 



fx^ds 



su PP (a-e4) = U r 2( s ) = r 2(0 



where the last equality holds since the sets ^(s) are decreasing as s increases, see 
also Propositions 11.41 and 11.51 

Thus in order to establish that (v\, v 2 ) is the minimizer of the energy functional 
(11.37)) under the conditions stated in Theorem 11.71 it suffices to prove that the 
variational conditions (|1.38[) and (|1.39[) are satisfied. 

The proof of (|1.38p and (|1.39[) will be carried out by integrating the variational 
conditions (|1.32[) and (|1.33p with respect to s from to £. The proof of Theorem 
2.3 of [6] contains a more general expression for the variational conditions (|1.32[) 
and (|1.33[) . which is valid for any complex number x, namely 



(5.4) 
(5.5) 



2 / log \x - y | df4(x) - J \og\x - y\d^ s 2 (x) -t = log 
2 J log \x - y\dfi 2 (x) - J \og\x- y\dfi°(x) = log 



Zi(x,s) 



z 2 (x,s) 
z 2 (x,s) 



z 3 (x,s) 



These conditions reduce to (|1.32[) and (| 1 .33[) whenever x E F-l(s) and x E r 2 (s), 
respectively. 

Proof of (jl.38p . If we multiply both sides of (|5.4|) by 1/^, integrate with respect to 

s from to £, and interchange the order of integration, then we obtain 

(5.6) 



2 J log\x- y\dv\{x) - J\og\x- y\dvl(x) -I = log 



z\(x, s) 



z 2 (x, s) 



ds, x E C, 



for some constant I E K, where v\ and v\ are the measures defined in (| 1 . 29[) and 
(|1.34p . respectively. 

Let x > 0. Since \zi(x, s)\ > \z2(x, s)\ for every s it follows from (15. 6p that 
(5.7) 



zi(x, s) 



z 2 (x, s) 



ds = -V{x). 



2 J log|x- y\dv{{x) - J log|x- y\dv^{x) -I < jj log 

If x G supp(^J) then x E T\(£) by (15. ip and therefore x E T\(s) for every s > £, 
since the sets are increasing. Thus |zi(x, s)\ — \z-z{x,s)\ for every s > and 
equality holds in (|5.7p for x E supp(i^). This completes the proof of (|1.38p . 
Proof of of (|1.39p . If we multiply both sides of equation (|5.5p by 1/^, integrate with 
respect to s from to £, and we interchange the order of integration we obtain 



log|x - y\dv 2 (x) - / log\x - y\dv((x) 



log 



zi{x, s) 
z 3 (x,s) 



ds. 
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Since \z2(x, s)\ > \zs(x, s)\ for every s > it follows that 

(5.8) 2 J log | a; - y\dv%(x) - J log \x - y\dvf(x) > 0, xeC. 

Equality holds in (|5.8[) if and only if \z2(x, s)\ = \za(x } s)\ for every s G (0,£], that 
is, if and only if 

x e P) r 2 (s) = r 2 (£) = supp(> - fz/J), 

where the first equality holds since the sets r 2 (s) are decreasing as s increases, and 
the last equality holds because of (|5.3p . This completes the proof of (|1.39[) . □ 



6. Proof of Theorem 11.81 

Before starting with the proof of Theorem ll.8l we establish the following lemma. 

Lemma 6.1. Let A s (z) be given by (11.27)) . and let Zi(x, s), Z2(x,s) andza(x,s) be 
the solutions of A s (z) = x, ordered as in (|1.15p . Then 



lim z\[x, s) = x — pt(l — t) — a(l — t) , 

s— >o+ 

,. — i / ^ 2at(l-t) 2 



I 6 - 1 ) *-o+ - 2a(l - t) + pt - ^pH 2 + 4ax 

lim s _1 Z3(x, s) = 



s ^°+ ' ' 2a(l -t)+pt- y/pH 2 + 4ax ' 

Proof. The lemma follows by a straightforward computation. □ 
For the proof of Theorem 11.81 we need to establish the two identities (|1.40p and 

Proof of (|1.40[) . Let x > 0. From Propositions 1 1 .41 and 1 1 . 51 it follows that there exist 
a unique s*(x) > so that for all s > 

x G I\(s) s > s*(x). 

Then log |zi(x, s)/z2(x, s)\ = for all s > s*(x), and so by (|1.35|) 

zi{x,s) 



(6.2) n, i = / log 

There is a special value 



z 2 (x,s) 



ds. 



x = (l- t)(a(l - t) +pt) = lim /3(s) = lim 7(5) 

that belongs to every Ti(s) for all s > 0. Then s*(xo) = and 
(6.3) V(x Q )=0. 
The derivative of (16.21) is 



( 6 .4) v (a0 = r x) ( 1 gfi^ii) l«) rfs . 

7 \zi(x, s) dx z 2 (x,s) dx J 
In order to handle this integral we introduce new variables 

, , ~, . zi(x,s) „. 22(2;, a) ~, \ z 3 (a;,s) 

(6.5) Zi(X,S) = , Z2{X,S) = , z 3 (x,s) = . 
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Since Zj(x,s) for j = 1,2,3 is a solution of A s (z) = x, it follows that Zj(x,s) for 
j = 1, 2, 3 is a solution of the equation B(z, s) = x where 

(6 6) B(z s) = ( sz + a ( 1 - t ) 2 )( z + t ( 1 - t )) 2 + Pt(l-t)(z + t(l-t) 
' z 2 z 

Taking partial derivatives with respect to s and x on both sides of B(zj(x, s), s) = 
x, and applying the chain rule, we obtain 

dB i~ i„ „\ „^^ 9zj(x,s) ( dB i 



— { Zj {x,s),s) J — t^— + —(zj(x,s),s) = 0, 

( 6 - ? ) ^ / Fiji \ a~ < y 

oB \ ozj(x, s) 



dz ^ ix > s) > s) )^^ = 1 
for j = 1, 2, 3. From (|6.6p . it is elementary to deduce that 

as Zj-(z, s) 

Combining this with (|6.7[) we obtain 

rfis^i 1 dz,(x,s) 1 dzj{x,s) 1 = 1 oo 

1 ' ^(a;,s) 9x (^(z,s)+t(l-t)) 2 9s ' J ' ' 

Using dnUI) in (JOJ) we get 

_ f<«> / 1 ^ l( x,s) _ 1 ^s) ' 

W 7o V (?i(a;,s)+<(l-t)) 2 9s (^(z, a) + t(l - i)) 2 9s J ' 

which can be written as 

(6.9) V'(x)=J ^(F(?i(a;,s))-F(z 2 (x,s)))rfs 

with 

1 



(6.10) F(x) = 

35 + t(l — l) 

From (|6.9p and the fundamental theorem of calculus, we have 

(6.11) V'(x) = F(zi(x,s*(x)) - F(Z 2 (x,s*(x)) - lim {F{z x {x, s)) - F(z 2 {x, s))). 

s— >0+ 

By definition s*(x) is the smallest value of s > for which x £ Ti(s). Then 
x = j(s*(x)) if xq < x and x — (3(s*(x)) if < x < Xo- We can observe from 
Figures [TTITOl that zi(j(s),s) — z 2 (7(s),s) and zi((3(s),s) — z 2 (/3(s), s). Therefore 

zi(x, s*(x)) = z 2 (x, s*(x)) 

and (|6.11[) reduces because of (16 . 1 0[) to 

(6.12) V'(x) = - lim f— — !— r-^ — ^ r ) . 

From (|6.5|) and Lemma 16. II we find that 

lim s) = oo and lim s) = ^ = ^ 

2a(l -t)+pt- v / P 2 i 2 + 4ax 

Then (|6.12p leads by straightforward computation to 

(6.13) V'(x)= ' 2 " 



t(l-t) t(y/ p H 2 +4ax-pt) 
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We obtain V(x) by integrating (|6.13|) with respect to x. Thus 
x 



V(x) = 



\/p 2 t 2 + 4ax , . — 

— p\og{yJp 2 t 2 + Aax - pt) + C. 



t{l-t) t 

The constant of integration C should be such that V(xq) = 0, see ()6.3|) . This leads 



to 



a= «a_j) +plog(2a(1 _ t)) 



and (jl.40p is proved. 

Proof of (|1.4ip . The measure a, introduced in (|1.36[) . has the density 

(6.i4) *p. = r*m d8i x<0 , 

ax J ax 
where /dx = for x ^ ^(s) and by (|1.30|) 

df4(x) 1 fz' 2 _(x,s) z' 2+ (x,s)\ _ 1 fz' 3+ (x,s) z' 2+ (x,s) 



dx 2m \Z2—(x, s) Z2+(x,s)J 2iti \Z3+(x, s) Z2+(x, s) 
for x G ^(s). The last equality comes from the fact that Z2-{x, s) — z 3+ (x, s) for 

x e r 2 (s). 

From Proposition 1 1 . 51 it follows that <i/if (&) /cfe = for every x G (— p 2 t 2 /4a, 0] 
and every s > so that 

^1 = forxG {-p 2 t 2 /4a,0}. 
dx 

Let x < p 2 t 2 /4a. Again from Proposition 11.51 (or from Proposition 11.41 in case 
p = 0) it follows that there is a unique s*(x) > so that 

x G ^(s) s < s*(x). 

Then d^, 2 (x)/dx = for s > s*(x) so that (]6 . 14[) is in fact a finite integral 

(6.15) ^ = j_r >} f i _g^£) i g^u, 

ctcc 2iri J \Z3+(X, s) ox Z2+{x,s) ox J 

The computation of (|6 . 1 5[) is similar to the computation of (|6 .4[) given above. 
We use the functions Zj(x,s) as in (|6.5[) and the function F as defined in (|6.10[) . 
From (|6.8p and (|6 . 1 5[) we then get 

Mx± r* {x) f i dz 3+ i dz 2+ \ 

dx J \(z 3+ +t(l-t)) 2 ds (I 2+ +t(l-t)) 2 ds J 

= F(z 3 +(x,s*(x))) - F(z2+(x,s*(x))) 

- lim (F(z 3+ (x, a)) - F(z 2 +(x, a))) 

(6.16) = - lim (F(z 3+ (x, a)) - F(z 2+ {x, a))) 

since z 2 +(x, s*(x)) = z 3+ (x,s*(x)). 

Since x < —p 2 t 2 /Aa, we obtain from (|6.5[) and Lemma HTT1 

lim S^fx, s) = lim 8~ 1 Z2(x, s) — } , 

•-0+ s^o+ 2a(l - t) + pt - iy/4a\x\ - P H 2 

lim z-i(x,s)= lim s _1 z-?(a;, s) — ; 

•^0+ v S -o+ 20(1 - 1) + pt + V 4a M - F 2 * 2 
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Using this and ()6.10|) in (|6.16|) we easily get 

„ d<j(x) ,. / 1 1 

2?ra r = - lim ' 



dx s->o+ \^2(a;, s) + t(l — t) z 3 (x, s) + t(l - t) 

. \f4a\x\ -p 2 t 2 

which completes the proof of identity (|1.4ip . □ 
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